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Abstract. We have performed a number of experiments with a Bose-Einstein 
condensate (BEG) in a one dimensional optical lattice. Making use of the small 
momentum spread of a BEG and standard atom optics techniques a high level 
of coherent control over an artificial solid state system is demonstrated. We are 
able to load the BEG into the lattice ground state with a very high efficiency by 
adiabatically turning on the optical lattice. We coherently transfer population 
between lattice states and observe their evolution. Methods are developed and 
used to perform band spectroscopy. We use these techniques to build a BEG 
accelerator and a novel, coherent, large-momentum-transfer beamsplitter. 
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An optical lattice is a practically perfect periodic potential for atoms, produced by 
the interference of two or more laser beams. A Bose-Einstein condensate (BEC) [Q, || 
is the ultimate coherent atom source, a collection of atoms, all in the same state, with 
an extremely narrow momentum spread. Combining a BEC with an optical lattice 
provides an opportunity for exploring a quantum system analogous to electrons in 
a solid state crystal but with unprecedented control over both the lattice and the 
particles. 

Periodic optical potentials have been widely used in atomic physics (For reviews 
see [S, 1^, pt). Recent experiments and theory have studied BECs in optical lattices 
§ Sliri^M HI HI' lil' II 111 lil' @ El' H' H' Hi H' H' H • Here we present 



a series of new experiments in which we precisely manipulate the BEC/lattice system 
and interpret the results in terms of a single-particle band structure theory. 

We describe ID lattice experiments with a sodium BEC in a regime where 
interactions between atoms are negligible. This paper is organised as follows: we begin 
with a description of the experimental arrangement (Sec. ^ and a brief summary of 
band theory (Sec. ||). We explore diabatic and adiabatic loading of a BEC into a 
lattice in Sec. ||. Using momentum state analysis we measure coherent superpositions 
between bands and determine that we can adiabatically load more than 99% of the 
atoms into a well-defined Bloch state. In Sec. ^ we manipulate the lattice potential to 
transfer population between Bloch states, mapping the band structure of the lattice 
eigenstates. We follow the work of |2^, |2|] in Sec. ^ in the deep lattice regime to 
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develop a BEC accelerator which imparts many photon recoils of momentum to the 
atoms while preserving the coherence of the condensate. With a little modification we 
turn the accelerator into a large momentum transfer beamsplitter, which is described 
in Sec. ^. 

1. Experimental setup 

We perform our experiments with a condensate of 3(1) x 10® sodium atoms |] in the 
35*1/2, F — fTT-F = —1 state. The sample has no discernible thermal component. 
The condensate is prepared as described in 1^ and is held in a magnetic time-orbiting 
potential (TOP) [^9[ |^ trap with trapping frequencies uj^ ~ \/2ujy = 2uJz — 27rx27 
Hz. Assuming a scattering length of a = 2.8 nm, the calculated Thomas-Fermi 
diameters |3^ are 47, 66, and 94 /im, respectively. The size and density of the BEC 
are such that the condensate populates on the order of a hundred lattice wells and 
interactions between atoms can be ignored on the time-scale of our experiments. 

After production of the condensate the optical lattice is turned on. It consists 
of two counter-propagating laser beams along the x-direction whose frequencies and 
amplitudes can be controlled independently by acousto-optic modulators (AOMs). 
The lattice beams are detuned about 60 GHz to the blue of the atomic resonance (the 
sodium D2 line: A — 2'K/k — 589 nm) so that the spontaneous emission rate (^ 100 
s~^) is negligible on the time scale of our experiments. The polarization is linear and 
parallel to the rotation axis of the TOP trap bias field (the y axis) . The condensate is 
located in the focus of the beams, which have a 1/e^ diameter of about 600 /im. The 
maximum power used in each beam is about 4 mW. 

The two counter-propagating laser beams form a standing wave, which acts on 
the atoms via the light-shift ||3l| to produce a sinusoidal potential: 

y(x,t) = y (H-cos(2fca; + 5t)). (1) 

Here 5 is the frequency difference between the two laser beams. This detuning gives 
the standing wave a velocity oiv = 5 /2k. A frequency difference of <5/(27r) = 100 kHz 
corresponds to a lattice velocity of 3 cm-s~^ (one photon recoil). 

We use time-of-flight analysis j^] to examine the momentum distributions that 
result from our experiments. At the conclusion of the experiment we turn the magnetic 
trap off ^ abruptly and the lattice is removed. The details of the turn-off of the lattice 
depend on the experiment and will be described later. Two milliseconds of free-flight 
allows the momentum components (which are discrete and separated by 21ik due to 
the periodicity of the lattice) to spatially separate. We then measure the number of 
atoms in each momentum state by absorption imaging with a CCD camera. 

We adjusted the laser power and detuning such that the height of the lattice 



potential Vb (measured as described in Sec. |3J|) was typically 14 photon recoil energies. 
En = fi'k'^ I2M. M is the atomic mass of sodium; Er/H — 25 kHz. The total power 
was constant to within 10%. 

II All uncertainties reported here are one standard deviation combined systematic and statistical 
uncertainties. 

^ It should be noted that, while the magnetic trap is on for the majority of the experiments, its 
presence does not affect our results. 
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2. Band structure in an optical lattice 

Band structure theory is derived assuming an infinite periodic potential. Our lattice 
period (w 0.3 fim) is small compared to the Thomas-Fermi diameter (47 /im) of the 
condensate in the lattice direction. Hence, the atomic wavefunction covers more than 
a hundred lattice sites and the system can be considered practically infinite. The 
corresponding momentum spread of the BEG is very small 0.01?ifc), so we treat 
the BEG as a plane wave. Also, atom-atom interactions are negligible: the time-scale 
associated with the peak chemical potential, h/fi = l/(27r x 800 Hz) — 280 /xs, is 
longer than the 20 /is time scale of our typical experiments. Thus, a single-particle 
picture is valid. 

The periodicity of the lattice leads to a band structure of the energy spectrum of 
the atom-lattice system. The eigenenergies En,q and eigenstates |n, q) (Bloch states) 
of the system in the rest frame of the lattice are labeled by the quasimomentum 
q and the band index n. The spatial periodicity of the eigenfunctions results in a 
momentum spectrum composed of peaks separated by the momentum corresponding 
to the reciprocal lattice vector, 2fik. The Bloch states jn, q) can thus be expanded in 
the discrete plane wave basis with momenta p — q -\- 2mhk : 



In order to calculate the coefficients an^q{m) and the eigenenergies -B^.g we solve the 
equation H\n, q) = En^q\n, q) where the Hamiltonian is H = P^/2i\/+(Vo/2) cos(2fca:). 

Figure ^ shows the eigenenergies versus quasimomentum for the first five bands 
{n — through 4) from q = —hk to Tik (the first Brillouin zone) for Vq = 14i?7?. Even 
in this relatively shallow potential the n — band is quite flat indicating that the 
tunneling rate is low. 

Figure Ijshows the plane wave decomposition, |a„.g(m)p, of the Bloch eigenstates 
for the bands n = 0, 1, 2, and 3 at q = and Vq — 14:Er. For example, the lattice 
ground state \n = 0, q — 0) consists mainly (65 %) of a momentum component and 
two smaller ±2hk momentum components (17 % each). 

The spatial wave function of a Bloch eigenstate in the lattice frame, tpn^qix) = 
{x\n,q), obeys Bloch's theorem: 



The quasimomentum characterizes the phase difference between neighbouring lattice 
sites. 

Optical lattices have the advantage that we can easily change their depth and 
velocity. In the absence of an external force (or equivalently an acceleration of the 
lattice), the quasimomentum is conserved (even if the depth changes) as long as the 
lattice periodicity is kept constant. (Recall that quasimomentum is only defined in the 
rest frame of the lattice.) For an optical lattice this is obvious since the only mechanism 
for changing momentum in the absence of spontaneous emission is redistribution of 
photons from one lattice beam to the other. Hence the momentum can only change 
by multiples of 2hk, which does not change the quasimomentum (see equation (||)). 



oo 




(2) 



m— — OO 



V'„,,(a: + A/2) = 6^9^/2'' ^„^,(a;) . 



(3) 
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Figure 1. The first five bands (even bands are solid, odd bands are dotted) of 
a lAEjf deep, ID, sinusoidal lattice across the first Brillouin zone. The dashed 
horizontal line indicates the top of the lattice. 



3. Loading the BEC into the lattice 

We now discuss methods for loading the BEC into the lattice with a selected 
quasimomentum q. We investigate two regimes: the sudden turn-on and the adiabatic 
turn-on. 

A BEC makes it easy to load a lattice state with a well-defined quasimomentum. 
Its small momentum spread (of much less than Uk) results in an similar spread of the 
quasimomentum q. For a condensate stationary relative to the lattice, the phase of 
the BEC is uniform across the lattice, hence q = 0. By contrast, when the condensate 
is moving with velocity v relative to the lattice, there is a linear phase gradient across 
the condensate and q = —Mv. 

3.1. Sudden loading of the lattice 

A BEC, taken to be a plane wave {\^{t = 0)) — \4>q)), suddenly loaded into a lattice 
can be described as a superposition of Bloch states \n,q) : 

oo 

\^{t = 0)) = J2\n,q){n,q\cl),) . (4) 

n=0 
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Figure 2. The lattice eigenstates can be expressed as a superposition of 
plane waves with fractional population \an,q(m)\^. Shown is the plane wave 
decomposition of the lattice eigenstates for the four lowest bands for q = and 
lattice height Vq = 14i?fl. 



From equation (g), {n,q\(j}q) = a* ,j(0). Therefore while the BEC wavepacket is held 
in the lattice it evolves in time according to 

oc 

I* W> = E {-»^«(<z)t/?i} \n, q) . (5) 

71=0 

We suddenly switch the lattice off after a time r. If we project the lattice state onto the 
plane wave (measurement) basis, we obtain the coefficients hq{m) of each \4>q+2mhk) 
in the lattice frame: 

^^<,q{^) an,q[m) CXp f -i -j^ T j . (6) 

The interference of the exponential factors produces oscillations in the populations of 
the plane wave components as a function of t. 

We observe these oscillations as follows: after producing the BEC we quickly 
(within 0.5 /is) turn on the lattice beams and leave them on for a variable length of 
time T. We then shut the lattice off abruptly (again within 0.5 /is) and measure, after 
a time of free-flight, the relative populations of the momentum components |^ 

+ This experimental set up is almost identical to the one published by Ovchinnikov et al. We 
discuss it here again for pedagogical reasons because Ovchinnikov et al. use a different model to 
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Figure 3. Coherent oscillations of the momentum decomposition of a BEC 
after suddenly switching on the lattice (g = 0, Vq = HEn). On the left are 
time of flight images showing the plane wave decomposition of the lattice state 
evolving as a function of time held in the lattice. The momentum components 
have been allowed to spatially separate. Only a single 4 /is cycle is shown. The 
right hand shows the populations of the 0, +2 and +4?ife momentum components 
(respectively: the filled squares, open squares and circles) over two cycles. The 
lines are a theoretical calculation of the populations in each component. 



Figure ^ shows the time evolution of a BEC held in the stationary lattice (q = 0). 
The populations vary almost purely sinusoidally, because only bands and 2 are 
significantly populated. Band 1 and all odd bands are not populated because the 
symmetry of those eigenstates are antisymmetric while the initial BEC wavefunction 
is essentially uniform (effectively symmetric) across a single well. The observed 
oscillation period of 250 kHz is the energy gap between band and 2 for a UEji 
deep lattice. We use this frequency to calibrate our lattice depth. 

We can also load the condensate into a lattice moving with a constant velocity 
relative to the BEC. In the rest frame of the lattice this corresponds to loading a 
specific non-zero q. Figure ^ shows the results of an experiment, similar to that shown 
in figure |^, for q = Ihk. The oscillations are not sinusoidal because more than two 
bands have significant population. Odd bands 1 and 3 are now populated in addition to 
bands and 2 because the selection rules discussed above are not valid ioi q ^ 0. The 
resulting beating signal consists of a superposition of sinusoids of different frequencies 
(see equation (^)) corresponding to the energy differences between pairs of bands. Our 
observations can be reproduced by the simple ID band model, as shown in the right 
hand side of figure ^. We have also performed the sudden loading experiment for other 
q between and Ihk and the results agree well with the mode]|^. 

We have measured the decay time of the beat signal by extending the lattice 
interaction time t. Figure ^ shows the decaying oscillations of the Ohk component for 
q ~ 0. Wc fit this data with an exponentially decaying sinusoid with a decay constant 
of 17 /iS. Notice that the data scatters widely at long times. This implies that the 

explain the underlying physics. Also we have extended that work from measuring at g = to 
arbitrary q and our results have significantly improved data quality. 

* Note that if we take q = Ihk and use a lattice shallow enough that only two bands are populated 
then we recreate the Bragg diffraction experiments described in M 
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Figure 4. Coherent oscillations of a BEC after suddenly switching on thfi 
lattice for q = ITik. This is the same as the experiment shown in figure H, 
only the quasimomentum has been changed. This time the oscillations are not 
sinusoidal. The right hand side plot shows the evolution of the —2hk component. 
The standard ID band model describes the data well, as shown by the full line. 



decay is due to two mechanisms: the inhomogeneity (« 10%) of our lattice beams over 
the condensate and laser power fluctuations from shot to shot (also about 10%). 

3.2. Adiabatic loading of the lattice 

In contrast to the preceding section, where we populated several bands by suddenly 
turning on the lattice, here we transfer the BEC into a single Bloch state \n,q). To 
achieve this we ramp up the lattice intensity adiabatically, i.e., such that p2| 

^AE^q,t)/h. (7) 

where AE is the energy difference between the ground state and the first excitable 
state \i). The left hand side is always less than dVb/dt. Therefore the condition in 
(0) can be satisfied at g = by choosing dVo/dt << WEj^/Ti since AE > 4Eu for 
any Vq (see figure]^). It becomes harder to maintain adiabaticity as q approaches the 
Brillouin zone boundary since there AE = when Vo = 0. The manifestation of the 
non-adiabaticity is Bragg scattering. For our experiments, at g = with a final lattice 
depth of lAEjj, an adiabatic ramp must be significantly longer than 5 /iS. 

Note that this adiabaticity criterion takes into account only the single particle 
motion. We are very non-adiabatic with respect to the time scale associated with 
atom-atom interaction energies. Recent experiments, which observe number squeezing 
and the Mott transition jl^, depend on being adiabatic on this longer time scale. 
If we can adiabatically ramp up the lattice intensity then we should also be able to 
adiabatically ramp down the intensity and return to the original BEC. This suggests 
an obvious test for adiabaticity where we ramp the intensity up and then down over 
equal time intervals. The degree of adiabaticity can be determined by measuring the 
number of atoms excited into momentum components other than hk. The results of 
this experiment are shown in figure |^. 

This method has two limitations: Firstly, the sensitivity is limited to the fraction 
of atoms that can be reliably detected, a few percent in our case. Secondly, interference 
between populated bands (as seen in the previous section) can lead to all the 
population being in Tik even though the process is far from adiabatic. Because 
of this, the signal oscillates and the adiabaticity of a single ramp time can only be 
determined by interpolating the envelope of the curve. 
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Figure 5. As the lattice is held on for longer times the amplitude of the 
oscillation signal decreases. This is due to two factors, inhomogeneity of the 
lattice beams and shot-to-shot variation in the lattice strength. Both of these 
effects lower the oscillation amplitude by averaging several oscillations which are 
slightly out of phase with each other. The inhomogeneity of the beams affects 
each data point and this decay is the decay of the envelope of the data points (this 
is indicated by the dashed lines which have been drawn to guide the eye). The 
shot-to-shot intensity fluctuations only have an effect when the data is averaged 
or a fit is made. This can be seen in the solid line which is a fit of a decaying 
sinusoid to the data points: the decay of the fit is faster than the decay of the 
envelope. Both these effects have to be taken into account when analysing the 
data. 

In order to directly measure the adiabaticity, i.e. the transfer efSciency into the 
lattice ground state, we perform experiments similar to those described in the previous 
section. Starting from a condensate at rest we ramp up the intensity of a stationary 
lattice||. We hold the BEC in the lattice for a time t, typically between to 10 fis, 
before suddenly switching off the light. We then study the oscillations in the plane 
wave decomposition of the lattice wavefunction as a function of r. Figure ^ shows 
an example of such oscillations for the Otik component at q = after a ramp time 
of 20 IIS. This beating signal has a very small amplitude compared to the beating 
signal of figure || for a lattice that was abruptly switched on. This indicates that 
most of the population is in band 0. If we had populated only the ground state, 
there would have been no beating at all. If we ignore the small oscillations, the 
measured momentum decomposition (60 % Ohk, 20 % +2tik, 20 % —2hk) is close to 

tt While we apply a linear voltage ramp, non-linearities in the response of the AOM smooth the ends 
(~10 % of the ramp time). 
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Figure 6. The energy eigenspectrum for two values of the quasimomentum, 
g = (left) and q = 0.3?iA: (right), as a function of the lattice height Vq. At 
Vb = the lattice eigenenergies are the free particle kinetic energies, p'^/2M and 
the eigenstates are the plane waves with momentum p = 2mhk + q (in the lattice 
frame) . 



the theoretical decomposition of the lattice ground state (Sec. |2|) . 

From the amplitude of the oscillations we can infer the population in the lattice 
ground state. To do this, we assume that only band and 2 are populated (since 
bands 1 and 3 cannot be excited at q = 0). Let p be the fraction of the population in 
band after loading into the lattice. Then the population in band 2 is 1 — p. For a 
lattice with a height of 14£'ij band has a calculated Ohk momentum component of 
65 % whereas band 2 has one of 34 % (see figure ||) . The fraction of population in the 
Ohk momentum component after evolving for a time r is 

P{t) = ^065^ + e'('"^+«) ^0.34(1 - p) ' (8) 

where to is the frequency difference between bands and 2 and is a phase that is 
dependent on the details of the ramp. The amplitude of the beating is 



AP = 2v/0.65p X 0.34(1 -p). (9) 

Equation (|^) does not take into account the decay of the beating amplitude due to 
dephasing. Before calculating p from the experimentally measured value for AP, we 
must first correct for the dephasing. We cannot directly use the decay rate obtained 
from figure ^ since, in this case, we use an intensity ramp rather than a constant 
intensity. To find the correct decay rate a simulation of the sudden turn on was 
performed with the effects of shot-to-shot noise and beam inhomogeneity included. 
The simulation was a direct solution of Schrodinger's equation with no mean-field 
effects. The parameters of the dephasing effects were then adjusted to match the 
observed decay and were then applied to a similar simulation which included the 
intensity ramp. 

Figure || shows how the oscillation amplitude decreases as a function of the lattice 
ramp time. The data are corrected for decay as described above. The largest correction 
is a factor of 1.8 increase in the oscillation amplitude. The data agree well with a 
calculation (full line) where the single-particle Hamiltonian corresponding to loading 
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Figure 7. Testing the adiabaticity of a linear intensity ramp with a double ramp. 
Perfect adiabaticity would leave all the atoms in the hk state. If the ramp is 
not fully adiabatic other momenta are populated. The observed oscillations in 
the ?i/jLcomponent are due to interference effects similar to the ones described in 
section w. The full line is a theoretical model obtained through direct integration 
of the Hamiltonian. It has been rescaled by about 10% in the time direction to 
match the data. 



a Vb = I'iEfj lattice is numerically solved. The lowest amplitude shown in figure ^ 
(the open circle, this data is shown in figure ^ corresponds to a ramp time of 20 
and a groundstate population of 99.7(1)% (the uncorrected result is 99.9%). 

Finally we point out that a linear ramp is not the most efficient way to populate 
the ground state. Numerical simulations ^ show that our actual ramp (smoothed by 
the AOM response) is more adiabatic than a linear ramp. These simulations agree 
well with our observed adiabaticity. Other ramp shapes could be made even more 
adiabatic. Since the band gap increases with the lattice intensity, the ramp can be 
accelerated as the intensity increases. 

As an alternative we find that we can dramatically improve the loading of the 
ground state by using interference effects. The single linear ramp can for example be 
broken up into two linear ramps with different slopes. Both ramps coherently drive 
population between the ground and excited states. The second ramp can then be used 



ft Carl Williams, private communication. 
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Figure 8. Residual oscillations in the plane wave decomposition of a lattice 
state loaded by slowly ramping up the lattice over 20/js. The beating is greatly 
reduced when we raise the lattice in this manner because the condensate has been 
loaded almost entirely into the lattice ground state. This data has been averaged, 
but not corrected for loss of amplitude due to inhomogeneity in the lattice beams 
or shot-to-shot variations in the lattice height. 
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Figure 9. The decrease of the beating signal with increasing ramp time for 
a final lattice height of lAEji. The data points are corrected to account for loss 
of coherence. The full line is from a band-theory calculation for a lattice with 
Vo = 14_Efl. The insert sketches how the lattice height is rajnped during the 
experiment. The open circle corresponds to the data in figure H. 
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to almost entirely cancel the effect of the first ramp|. Details of these studies will be 
subject of future work. 

4. Coherent Transfer between Lattice Bands and Band-spectroscopy 

Starting from a condensate adiabatically loaded (20 fis ramp) into the ground state 
of the lattice, we will now study coherent population transfer to higher bands. 

4.I. Transfer between bands and 1 

To transfer population from band to 1 at g = (see figure |l^) we modulate the 
phase of one of the lattice beams and thus the position of the wells, i.e. shaking the 
lattice. The modulation is performed using an AOM and has an amplitude of 7r/6. 
We apply the modulation for four cycles over approximately 30 /is. 

Shaking the lattice is an odd-parity excitation and therefore efficiently couples 
bands of opposite parity, such as bands and 1. The phase modulation can also 
be viewed as putting sidebands on the lattice beam. In this picture the transition 
between bands is a Raman transition and the ability to make odd parity excitations 
is related to the phase difference of tt between the sidebands. 

A convenient way to analyze the population of excited bands is to adiabatically 
ramp down the lattice after the shaking. Population in band 1 at q = is all 
transferred, after the ramp-down, into the ±2hk components. The frequency of the 
to 1 transition can then be found by maximising the population in these momentum 
components. 

Once the resonance is found we can examine the population transfer via 
interference experiments. If we switch off the lattice suddenly we observe oscillations 
as a function of the holding time (see figure |l^ B and C). Since bands and 1 have 
opposite parity the oscillations will tend to be between momentum components with 
opposite momenta (in contrast to the beating of bands and 2 illustrated in figure 
Fitting an exponentially decaying sinusoid to the oscillating momentum components 
yields an oscillation frequency of 140 kHz for the ±2?ifc components consistent with 
the calculated frequency of 147 kHz for the llEji lattice used for this experiment. The 
oscillations show that we have coherently populated band 1. The small oscillation of 
the Tik momentum peak is due to population in band 2 (the oscillation frequency 
corresponds to the energy gap between bands and 2) that has been excited by 
second-order processes. 

4-2. Transfer between bands and 2 

We can transfer population between bands and 2, without exciting band 1, by 
amplitude-modulation of the lattice, an even-parity excitation. This can be viewed as 
a Raman process as well; in this case the sidebands are in phase. 

Although band 1 is relatively flat, band 2 has significant curvature (see Figs. |l| and 
|ll| A) which can be examined by varying q and measuring the resonance frequency. 
The resonance is found by varying the modulation frequency to maximise the loss 
from band 0. 

f Our numerical calculations show that using two optimized linear ramps one can achieve a ground 
state population greater than 99.98 % for a total ramp time of 20 /^s and Vq = liEji. This value for 
the ground state population is even stable to lattice intensity fluctuations of 10%. 
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Figure 10. Starting from the lattice ground state (q = 0, n = 0) population 
can be coherently transferred to the band 1 (at g = 0) by shaJfing the lattice for 
30 fis. This is illustrated in A. In B and C the oscillations in the plane wave 
decomposition of the lattice state are shown after transferring some population 
to band 1. Oscillations are strongest between the ±2%k components. Fitting 
exponentially decaying sinusoids (C) to the oscillating momentum components 
yields an oscillation frequency of 140 kHz for the ±2%k components. 

To produce a non-zero q within the first Brillouin zone we start from an 
adiabatically loaded lattice ground state {q = 0) and then suddenly apply a constant 
acceleration a. = 300 m-s~^ to the lattice. For |Ma(l, 0|a:;|0, 0) | ^ AE (corresponding 
to a < 3.9 X 10^ m • in our case) this sudden change does not significantly 
excite higher bands. In the lattice frame the inertial force changes q according to 
dq/dt = Ma. Wc adjust the duration of the acceleration to obtain the desired q. 
(We could have adiabatically loaded the condensate directly into a moving lattice to 
produce a non-zero q, but achieving adiabaticity is difficult near q = Ihk.) 

We suddenly stop accelerating, holding q constant. The lattice height is then 
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Figure 11. The scheme for transferring population into band 2 at arbitrary q. 
We start from the ground state at a quasimomentum of 0. We then accelerate the 
lattice to change the quasimomentum from to the desired value of q. Population 
is transferred into band 2 by modulating the lattice height over a 20 to 30 fis 
period. In order to measure the population in the upper band the acceleration is 
reversed to bring the quasimomentum back to and the lattice is adiabatically 
released to map the bands into individual momentum components. 



modulated by 15% (controlled with an AOM) for 6 to 7 oscillations. Using this 
relatively short and strong agitation to transfer the population broadens the transition 
and makes it easier to find the resonance line. In order to measure how much 
population has been transferred to band 2, the lattice is decelerated to zero velocity 
(g — 0) and then adiabatically ramped down to map the lattice bands into plane 
waves (see figure ^). Figure |l^ shows the resonance curve for q = Ihk. The resonance 
frequency at about 320 kHz equals the expected band gap between band and 2 
at q = Ihk. The width of the excitation is determined by the bandwidth of the 
excitation (« 37 kHz). In the insert we plot the measured resonance frequencies for 
the quasimomenta g = 0, 0.25, 0.5, 0.75 and Ihk together with the calculated band 
structure for Vq = UEj^. The data agrees well with the band model. 

5. An accelerator for a BEC 

An accelerated lattice can be used to reach velocities corresponding to momenta well 
beyond the first Brillouin zone (Similar work with thermal atoms has been described 
in 11^]). Using the method described above we adiabatically load the BEC into a 
deep lattice with a 25 ramp. We then accelerate the lattice at 1600 m-s"^. 
This sweeps the quasimomentum through a Brillouin zone every 35 /iS. 

Figure |l^ shows the change of the mean momentum and the momentum 
distribution of the accelerated BEC in the laboratory frame measured by suddenly 
switching off the lattice after a given period of acceleration. We see a linear increase 
of the mean momentum (corresponding to the lattice velocity) . Note that despite the 
continuous change of the mean momentum, individual atomic momenta (measured in 
the lab frame) are always multiples of 2hk. After passing through each Brillouin zone 
momentum decompositions are similar and recurring, shifted only by multiples of 2hk. 
This recurrence is due to the recurrence of q as it scans the Brillouin zone from ~~hk 
to hk and then undergoes Bragg scattering back to —hk. 

Despite the similarity of our experiments to those which exhibit Bloch oscillations 
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Figure 12. Experimental data showing population transfer from band to band 
2 at g =J£k as a function of excitation frequency. The scheme used is outlined 
in figure hll. The fit is a Lorentzian. The resonance frequency corresponds to the 
band gap between band and 2 at g = Ihk. Similar measurements at g = 0, 0.25, 
0.5 and 0.75^fc allow us to map out the curved band structure (see insert). 



in the group velocity [g6|, |2^, |l2| the osciUations in our case are too smaU to be 
observed. The group velocity of the atoms in the lattice frame (given by Vg = dE/dq) 
is, in our case, at most 7.3x10^^ m-s^^, which is negligible compared to the lattice 
velocity. (The maximum lattice velocity shown is 2.1x10^^ m-s""'^.) With increasing 
potential the Bloch oscillations in the group velocity are suppressed and we are left 
only with the oscillations in the momentum decomposition. 

To obtain an accelerated BEC with a single momentum we turn the lattice off 
adiabatically. The BEC has a final momentum of 2mhk where m is the number of 
Brillouin zone boundaries crossed during the acceleration (i.e. the final lattice velocity 
is between {2m — l)hk/M and (2m + \)hk/M.). We have accelerated a BEC up to 
momenta of IQfik (see figure |l^ B). The final velocity is quite insensitive to power 
and frequency fluctuations in the lattice beams and to the specifics of the lattice 
acceleration. The limitations on our final velocity are purely technical: the range of 
the frequency sweep and the field of view of our camera. 

While there is no apparent limit on the final velocity there are limits on the 
acceleration. In our case the sudden initiation of the acceleration presents a limit as 
described above. If we adiabatically applied the acceleration to circumvent this limit 
there will always be the ultimate limit of "spilling" : when the acceleration is so large 
that there are no bound states. 




20 40 60 80 100 120 

Time (|lIs) 

Figure 13. Tlic momentum decomposition and mean momentum of a 
condensate in the lowest band after a period of constant a<;celeration. The 
acceleration sweeps the atoms through a Brillouin zone in 35;us. The reference 
frame is the stationary laboratory. 
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Figure 14. Accelerator for cold atoms. First the BEC is adiabatically loaded 
into the lattice ground state {q =0) and the lattice is subsequently accelerated 
(A). After sweeping through an integer number m of Brillouin zones, the lattice 
is ramped down adiabatically, which results in a single BEC wavepacket moving 
at 2mhk. (B) shows accelerated BEC wavepackets with final momenta of up to 
lOhk. 



6. A large-momentum-transfer beamsplitter 

We now describe liow we use ttie accelerator in order to build a coherent, large- 
momentuni-transfer beamsplitter (LMT-beamsplitter). This is useful for applications 
that require atoms with large, coherent separations (in either position or momentum). 
A LMT-beamsplitter can be built using high-order Bragg diffraction j|] , but is limited 
to relatively low momenta because the transition probability decreases rapidly with 
increasing order. Our approach is not limited in this way. It is similar in spirit to 



coherent beamsplitters based on adiabatic transfer |3^, |3J, 35 . 

Our scheme is outlined in figure |l^: We begin with a second-order Bragg pulse 
P that produces a coherent superposition of the two momentum states and -Ahk. 
We use a SEu high lattice pulsed on for 6.4 iis at a velocity of —2hk/M to perform the 
splitting. We then adiabatically turn on a 14i?fl deep lattice at a velocity of 0.3hk/M. 
The lattice is accelerated to a velocity of approximately 2mhk/M carrying the Ohk 
momentum component along with it. By contrast, the -Ahk momentum component 
remains stationary. 

The selective acceleration is achieved by choosing a lattice velocity such that the 
atoms with momentum Ohk in the lab frame are in band 0, while the atoms with -Ahk 
are in band 4. Both have a quasimomentum of —0.3hk. We chose this q to ensure 
there are no degeneracies between bands at Vq = (see right hand side of figure |^) . 

As we accelerate the lattice, the atoms in band are accelerated with the lattice, 
just as described in the previous section. Crossing the Brillouin zone boundary does 
not change the band. By contrast the atoms in band 4 respond as nearly free particles 
(since band 4 is unbound) : when q reaches the avoided crossing between bands 4 and 5 
at the zone boundary (see figure |^) the atoms pass through to band 5 diabatically. At 
subsequent level crossings (such as for bands 5 and 6 at g = 0) this process continues. 
In the lab frame this corresponds to the atoms being unaccelerated. 

We can apply the Landau-Zener criterion for diabatic passage through the avoided 
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Figure 15. A coherent, large-momentum, transfer beamsplitter. After splitting 
the condensate with a second-order Bragg pulse, a lattice moving at 0.3hk/M 
is adiabatically turned on over 20 fis, such that the Ohk and —4hk components 
are transferred into bands and 4, respectively. The lattice is then accelerated 
over 200 fis. The component in band accelerates with the lattice, while the 
component in band 2 stays at rest. After sweeping through m Brillouin zone 
boundaries, the lattice is adiabatically turned off, resulting in two wavepackets, 
one with momentum -4hk and one with a momentum of 2mhk. In our experiments 
we swept through up to 4 Brillouin zone boundaries in about 200 fis. 

crossings. The probability for a passage between bands (n — 1) and n is [|7j : 

=cxp|-7r-^V (10) 

where w is the band gap frequency. This expression is derived in the Umit where the 
anti-crossings are small. For the passage from bands 4 to 5 the gap is 2.2 kHz, giving a 
diabatic probability of 99.8% for our acceleration of a = 1200 m-s~^. (If we had used 
a first-order Bragg pulse to initially split the condensate into and 2fik components, 
then only 72% of the population would have passed through the first crossing.) 

We end the experiment by adiabatically ramping down the intensity of the moving 
lattice. We get two wave packets, one at 2mhk (recall that ni is the number of Brillouin 
zone boundaries crossed), the other at rest. We have demonstrated splittings of up to 
12hk. 

The LMT-beamsplitter can be used to construct an interferometer. We have 
performed a proof-of-principle experiment in which we have built a Mach-Zehnder 
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interferometer similar to that described in |3^, In our case the initial Bragg 

pulse beamsplitter was replaced with the LMT-beamsplitter and the final Bragg pulse 
beamsplitter with a time-reversed version of the LMT-beamsplitter. In addition the 
Bragg TT-pulse of the conventional Mach-Zehnder interferometer is replaced with a 
three stage LMT-7r-pulse consisting of a deceleration stage, a conventional Bragg tt- 
pulse and a final acceleration stage. The resulting interference pattern had a contrast 
of 50%, which was limited by transverse motion of the BEC during the experiment. 

7. Conclusions 

In conclusion, we have performed a number of experiments with a BEC in static, 
moving, and accelerating ID optical lattices and have interpreted the results using 
band structure. We have studied adiabatic loading into the ground state and coherent 
transfer to other lattice states by modulating the lattice. These techniques have been 
used to perform band spectroscopy. We have also built an coherent accelerator and a 
multiphoton beamsplitter. 

During preparation of this paper we learned of several other optical lattice 
experiments with a BEC [|l| |l|, |l|, |l|, |l^, |l|] . 

Acknowledgments 

We thank Keith Burnett, Mark Edwards, Lu Deng and Ed Hagley for fruitful 
discussions. This work was supported in part by the US Office of Naval Research, 
NASA and ARDA/NSA. J. H. D. and H. H. acknowledge support from the Alexander 
von Humboldt foundation. A. B. acknowledge partial support from DGA. 

References 

[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell. Science, 
269:198, 1995. 

[2] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Drutcn, D. S. Durfce, D. M. Kurn, and 

W. Ketterle. Phys. Rev. Lett, 75:3969, 1995. 
[3] J.P. Dowling and J. Gea-Banacloche. Adv. At. Mol. Opt. Phys., 37:1, 1996. 
[4] R. Grimm and M. Weidemiiller. Optical dipole traps for neutral atoms. Adv. At. Mol. Opt. 

Phys., 42:95, 2000. 

[5] S. Bernet, R. Abfalterer, C. Keller, M. K. Oberthaler, J. Schmiedmayer, and A. Zeilingor. Phys. 

Rev. A, 62:023606, 2000. 
[6] B.P. Anderson and M. A. Kasevich. Science, 282:1686, 1998. 

[7] Yu. B. Ovchinnikov, J. H. Miiller, M. R. Doery, E. J. D. Vredenbregt, K. Helmerson, S. L. 

Rolston, and W. D. Phillips. Phys. Rev. Lett., 83:284, 1999. 
[8] M. Kozuma, L. Deng, E. W. Hagley, J. Wen, R. Lutwak, K. Helmerson, S. L. Rolston, and 

W. D. Phillips. Phys. Rev. Lett., 82:871, 1999. 
[9] J. Stenger, S. Inouye, A. P. Chikkatur, D. M. Stamper-Kurn, D. E. Pritchard, and W. Ketterle. 

Phys. Rev. Lett., 82:4569, 1999. 
[10] L. Deng, E. W. Hagley, J. Denschlag, J. E. Simsarian, M. Edwards, C. W. Clark, K. Helmerson, 

S. L. Rolston, and W. D. PhiUips. Phys. Rev. Lett., 83:5407, 1999. 
[11] W. K. Hensinger, H. HafTner, A. Browaeys, N. R. Heckenberg, K. Helmerson, C. McKenzie, G. J. 

Milburn, W. D. Phillips, S. L. Rolston, H. Rubinsztein-Dunlop, and B. Upcroft. Nature, 

412:52, 2001. 

[12] O. Morsch, J. H. Miiller, M. Cristiani, D. Ciampini, and E. Arimondo E. Phys. Rev. Lett., 
87:140402, 2001. 

[13] F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. Smerzi, and 
M. Inguscio. Science, 293:843, 2001. 



A Bose-Einstein Condensate in an optical lattice 



20 



[14] P. Pedri, L. Pitaevskii, S. Stringari, C. Fbrt, S. Burger, F. S. Cataliotti, P. Maddaloni, F. Minardi, 

and M. Inguscio. Phys. Rev. Lett, 87:220401, 2001. 
[15] S. Burger, F. S. Cataliotti, C. Fort, F. Minardi, M. Inguscio, M. L. Chiofalo, and M. P. Tosi. 

Phys. Rev. Lett, 86:4447, 2001. 
[16] M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and T. Esslinger. Phys. Rev. Lett, 87:160405, 

2001. 

[17] C. Orzcl, A. K. Tuchman, M. L. Fcnsclau, M. Yasuda, and M. A. Kascvich. Science, 291:2386, 
2001. 

[18] M. Greiner, O. Mandel, T. Esslinger, T. W. Hanscli, and I. Bloch. Nature, 415:39, 2002. 
[19] K. Bcrg-S0rcnscn and K. M0lmer. Phys. Rev. A, 58:1480, 1998. 

[20] D. Jalcsch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller. Phys. Rev. Lett, 81:3108, 
1998. 

[21] D. Clioi and Q. Niu. Phys. Rev. Lett, 82:2022, 1999. 
[22] M. Holtliaus. J. Opt. B, 2:589, 2000. 

[23] M. L. Chiofalo, M. Polini, and M. P. Tosi. Eur. Phys. J. D, 11:371, 2000. 

[24] S. Potting, M. Cramer, C. H. Schwalb, H. Pu, and P. Meystre. Phys. Rev. A, 64:023604, 2001. 
[25] J. C. Bronski, L. D. Carr, B. Deconinck, and J. N. Kutz. Phys. Rev. Lett, 86:1402, 2001. 
[26] M. Ben Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon. Phys. Rev. Lett., 76:4508, 1996. 
[27] E. Peik, M. Ben Dahan, I. Bouchoule, Y. Castin, and C. Salomon. Phys. Rev. A, 55:2989, 1997. 
[28] S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Qian Niu, and M. G. Raizen. Phys. Rev. 
Lett, 76:4512, 1996. 

[29] W. Petrich, M. H. Anderson, J. R.Ensher, and E. A. Cornell. Phys. Rev. Lett, 74:3352, 1995. 
[30] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari. Rev. Mod. Phys., 71:463, 1999. 
[31] H. J. Metcalf and P. van der Straten. Laser Cooling and Trapping. Springer, 1999. 
[32] L. SchifT. Quantum Mechanics. McGraw Hill, third edition, 1968. 
[33] J. Lawall and M. Prentiss. Phys. Rev. Lett, 72:993, 1994. 

[34] L. S. Goldner, C. Gerz, R. J. C. Sprceuw, S. L. Rolston, C. I. Westbrook, W. D. Phillips, 

P. Martc, and P. Zoller. Phys. Rev. Lett., page 997, 1994. 
[35] M. Weitz, B. C. Young, and S. Clm. Phys. Rev. A, 50:2438, 1994. 

[36] P. J. Martin, B. G. Oldaker, A. H. Miklich, and D. E. Pritchard. Phys. Rev. Lett., 60:515, 1988. 
[37] G. Zcncr. /Voc. R. Soc. London A, 137:696, 1932. 

[38] J. Denschlag, J. E. Simsarian, D. L. Feder, C. W. Clark, L. A. Collins, J. CubizoUes, L. Deng, 
E. W. Hagley, K. Helmerson, W. P. Reinhardt, S. L. Rolston, B. I. Schneider, and W. D. 
Phillips. Science, 287:97, 1999. 

[39] Y. Torii, Y. Suzuki, M. Kozuma, T. Sugiura, T. Kuga, L. Deng, and E. W. Hagley. Phys. Rev. 
A, 61:041602, 2000. 



